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Abstract. We investigate the energy dependencies of X-ray quasi-periodic oscillations in black 
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hole X-ray binaries. We analyze RXTE data on both the low- and high-frequency QPO. We 
construct the low-/ QPO energy spectra, and demonstrate that they do not contain the thermal 
disk component, even though the latter is present in the time averaged spectra. The disk thus 
^ does not seem to participate in the oscillations. Moreover the QPO spectra are harder than the 
^ O time averaged spectra when the latter are soft, which can be modeled as a result of modulations 
occurring in the hot plasma. The QPO spectra are softer than the time averaged spectra when 
the latter are hard. The absence of the disk component in the QPO spectra is true also for the 
high-frequency (hecto-Hz) QPO observed in black hole binaries. We compute the QPO spectra 
I I expected from the model of disk resonances. 
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, 1. Introduction 

\^ Quasi-periodic components of X-ray variability are often seen in power density spectra 

' (PDS) from accreting compact objects, signifying the existence of "clocks" operating on 

I a specific time scale, among the broad-band variability. The nature of these periodic 

rS i' processes is unknown, but they are usually associated with various modes of oscillations 

^ ' of standard accretion disks present in those systems. 

^ Energy dependencies of the QPO can be an important clue to the origin of these 

, features, similarly to analyzes of time averaged energy spectra, which provide informa- 

Ci • tion about physical processes of generation of X-rays. It was noted quite early on into 

^ , the study of QPO that their r.m.s. variability seems to increase with energy (review in 
van der Klis (2006)] ). Here we report on more detailed investigations of QPO energy de- 
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pendencies in black hole X-ray binaries. Firstly, we construct observed low-/ (1-10 Hz) 
QPO energy spectra in the 2-20 keV range using RXTE data. Secondly, we construct 
models of spectral variability of inverse-Compton emission, producing a given type of 
QPO spectra. We also analyze data on high-/ (hecto-Hz) QPO and construct models of 



energy spectra of these QPO in the model of disk resonances ( Abramowicz & Kluzniak (2001) | 



2. Low-frequency QPO 

2.1. Observations 

We analyzed observations of four sources: GRS 1915-1-105, XTE 1550-564, 4U1630-47 
and XTE 1859-f 226 performed by RXTE in observing modes providing good energy 
and timing resolution. Details of data reduction and analysis procedures are given in 
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Figure 1. Comparison of time averaged spectra and the QPO spectra (upper panels. Lower 
panels show ratios of the QPO spectra to various components from the time averaged spectra. 



Sobolewska & Zycki (2006) Firstly, we construct time-averaged energy spectra in the 
3-200 keV energy range using PCA and HEXTE data. These are fit with a model con- 
taining two or three continuum components (disk black body, comptonization) and the 
reprocessed component (Fe Ka line and the Compton reflected continuum). Then we 
construct the QPO energy spectra: we compute PDS in each energy channel, then fit 
it with a model consisting of a broad band continuum (broken power law or a number 
of Lorentzians) and one or two narrow Lorentzian peaks representing the QPO, possi- 
bly with a harmonic. The QPO Lorentzians are then integrated over frequency. This, 
expressed as variance (rather than the usual r.m.s./mean) is the QPO energy spectrum. 
The meaning of such a function is that it would correspond to energy spectrum of a 
variable component, if it was a separate spectral component which would be responsible 
to the QPO. 

The results are presented in Fig [1] We show both the time averaged spectra and the 
QPO spectra. We also show the ratios of QPO spectra to the time averaged models. 
The main result is the lack of the soft disk component in the QPO spectra. Presence 
of this is excluded in all analyzed spectra at high confidence level. This can be seen for 
example in the spectra of XTE 1859-1-226, where the QPO spectra are much closer to the 
Comptonized component only than to total time averaged spectra. Secondly, we analyzed 
the slope of the QPO spectra relative to the time averaged spectra. In soft spectral state 
the QPO spectra are clearly harder than the time averaged spectra (Fig[T^), while the 
opposite is true for the hard state (Fig[TJ;). During the two observations of 4U1630-47 
the time averaged spectra are intermediate in slope between values typical for the soft 
and hard states (F 2). We find that the QPO spectra are similar to the time averaged 
spectra in these two observations. 

2.2. Models 

In an attempt to understand the results presented in the previous Section, we performed 
simulations of spectral variability of inverse-Compton component. Such spectra are de- 
scribed by three main parameters: the plasma heating rate, /h, the flux of cooling soft 
photons. Is, and the temperature of the input soft photons, Tq. Irrespectively on the 
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Figure 2. Model predictions for tlie QPO energy spectra (red circles) in various scenarios for 
modulations: a) heating, b) cooling, c) covering factor of cold matter. Green squares are the 
broad band noise spectra, blue stars: QPO harmonic, solid line shows time averaged spectra. 



physical mechanism of QPO production, it is one of the three parameters that has to 
be modulated if the observed QPO corresponds to modulation of intrinsic luminosity 
(other possibility is a periodic modulation of geometry affecting, e.g., relativistic effect 
influencing the observed luminosity; see below). 

The slope of the time averaged Comptonized spectrum depends (almost) only on the 
ratio Ih/ls, however the mode of spectral varia bility depends on which para meter is modu- 
lated. Details of this investigation are given in Zycki & Sobolewska (2005) To summarize 



it: modulation of the heating rate produces variability spectrum with largest amplitude of 
variability at high energies, which means that the variability spectrum is harder than the 
time averaged spectrum. Modulation of the cooling rate or the soft photons temperature 
produces variability spectrum softer than the time averaged spectrum. We considered 
also modulation of the amplitude of the reprocessed component, which was assumed to 
relate the cooling and heating rate (see, e.g., Poutanen & Fabian (1999) In this case we 
also find that the QPO are softer than the time averaged spectra, but, additionally, a 
very strong Fe Ka line is seen in the QPO spectrum. 

Summarizing, properties of low-/ QPO spectra in the soft state strongly suggest that 
these QPO are generated as a result of modulations in the hot plasma rather than 
in the standard cold accretion disk. This is contrary to most theoretical ideas where 
oscillations of the cold disk are thought to be the primary driver. However, the physical 
mechanism does seem to require both accreting plasma phases to work, since these QPO 
are associated with changes of the geometry of the system. It is also interesting to note, 
that long-time scale variability of Cyg X-1 in the soft state also seems to be driven by 
instabilities in the hot Comptonizing plasma (Zdziarski et al. (2002)). 



3. High-frequency QPO 

Data on high-frequency QPO in black hole X-ray binaries are of much poorer qual- 
ity than data for the low-/ QPO. It is nevertheless clear that they share at least one 
characteristics with their low-/ counterparts, namely the soft disk component is ab- 
sent in the QPO spectra. More detailed energy dependencies (i.e., the spectral slope) 
are not possible to determine with current data. The fact that pairs of QPO with 3:2 
frequency ratio are sometimes observed, stimulated development of models invoking non- 
linear resonances between various modes of oscillations of accretion disks. In particular. 



Abramowicz & Kluzniak (2001) postulate that the relevant modes are the vertical and 
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Figure 3. QPO spectra from the model of epicyclic oscillations, with relativistic effects mod- 
ulating the X-ray flux, for a — 0.998 and inclination /lobs ~ 0.3. Left: geometry of accretion 
disk with an active corona produces strong disk component in the QPO spectra. Right: hot 
oscillating torus produces QPO spectrum resembling a comptonization component. 

radial epicyclic motions. Importantly, they suggest also that the X-ray flux modulation 
occurs due to changes of geometry and, consequently, changes in the strength of the 
relativistic efl^ects (Doppler and gravitational shifts; light bending), with the primary 
emission being constant. This is thus the only current model making definite predictions 
for the mor e detailed properties of the Q PO, for example, their energy spectra. We used 



the code of Zycki & Niedzwiecki (2005) for photon transfer in Kerr metric and modified 
it to implement the source motion corresponding to disk epicyclic oscillations (details are 
in Zycki et al., in preparation). We have performed computations of the relativistic effects 
expected from such disks at radii corresponding to 3:2 resonance between the oscillation 
frequencies, and we have computed energy spectra of QPO expected in this model. They 
are presented in Fig. [3] in two variants of geometry: (1) accretion disk extending to the 
last stable orbit with an active X-ray corona and (2) a truncated accretion disk with inner 
hot Comptonizing flow. In the first scenario the QPO spectrum contains a strong disk 
component, since the disk emission is modified by the relativistic effects. In the second 
scenario the disk is, by assumption, not present at the QPO radius and, as a consequence, 
only the Comptonized component is modulated. The Fe Ka line is also assumed to be 
produced at the QPO radius (in the accretion disk with a corona geometry) and it does 
appear in the QPO spectrum, at positions corresponding to maximum redshift expected 
from the oscillatory motion. 
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